Few-particle quantum dynamics — comparing 
Nonequilibrium Green's functions with the 
generalized Kadanoff-Baym ansatz to density 
operator theory 

S Hermanns 1 , K Balzer 2 and M Bonitz 1 

1 Institut fur Theoretische Physik und Astrophysik, Universitat Kiel, D-24098 Kiel, Germany 

2 University of Hamburg, Max Planck Research Department for Structural Dynamics Building 
99 (CFEL), Lumper Chaussee 149, D-22761 Hamburg, Germany 

E-mail: hermanns@theo-physik.uni-kiel .de 

Abstract. The nonequilibrium description of quantum systems requires, for more than two 
or three particles, the use of a reduced description to be numerically tractable. Two possible 
approaches are based on either reduced density matrices or nonequilibrium Green's functions 
(NEGF). Both concepts are formulated in terms of hierarchies of coupled equations — the 
Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy for the reduced density operators 
and the Martin-Schwinger-hierarchy (MS) for the Green's functions, respectively. In both 
cases, similar approximations are introduced to decouple the hierarchy, yet still many questions 
regarding the correspondence of both approaches remain open. 

Here we analyze this correspondence by studying the generalized Kadanoff-Baym ansatz 
(GKBA) that reduces the NEGF to a single-time theory. Starting from the BBGKY-hierarchy we 
present the approximations that are necessary to recover the GKBA result both, with Hartree- 
Fock propagators (HF-GKBA) and propagators in second Born approximation. To test the 
quality of the HF-GKBA, we study the dynamics of a 4-electron Hubbard nanocluster starting 
from a strong nonequilibrium initial state and compare to exact results and the Wang-Cassing 
approximation to the BBGKY hierarchy presented recently by Akbari et al. [T]. 



1. Introduction 

The ab initio time-dependent description of quantum many-body systems has been a major focus 
in the physics and chemistry communities since the invention of quantum mechanics. Despite 
many efforts, up to the present time, the exact (analytic or numerical) solution of the underlying 
equation of motion for the wavefunction \^?) s (t) of the system — the Schrodinger equation — is, 
in general, restricted to only very few particles, due to the exponentially growing complexity 
with the particle number, e.g. [2J. To overcome this limitation a great variety of methodologies 
have been developed to approximately describe systems of larger particle numbers. This includes 
time-dependent density functional theory |3J , time-dependent Hartree-Fock [H El El E] , multi- 
configuration time-dependent Hartree-Fock [9], time-dependent coupled cluster theory |1(J] . 
M0ller-Plesset many-body perturbation theory [TT] , density matrix renormalization group based 
approaches [T2] or dynamical mean field theory, e.g. |13j . 



In this article, we focus on two different closely related powerful methods: non-equilibrium 
Green's functions (NEGF) [JJ] and reduced density operators (DO), e.g. [TdI [JjJl El US], 
respectively. Both methods involve the solution of a coupled hierarchy — the MS and BBGKY 
hierarchy, respectively — of reduced quantities, where the solution of the full hierarchy is equivalent 
to the solution of the iV-particle Schrodinger equation or von Neumann equation, in the case of 
pure or mixed states, respectively. To reduce the complexity, both methods aim at simplifications 
via a physically motivated truncation of the hierarchy. Since both approaches are selfcontained and 
independent, comparisons of the two are of great interest for applications. NEGF are commonly 
regarded as more accurate, however, nonequilibrium solutions of the two-time equations of 
motion — the Keldysh-Kadanoff-Baym equations (KBE) — are computationally very expensive. 
While there has been remarkable progress during the last decade for homogeneous |19[ I2L)| 
and inhomogeneous systems, e.g. [251 E3] , the two-time structure puts strong limits on the 
achievable propagation time. In contrast, the single-time density operator approach does not 
suffer this problem. Furthermore, recent solutions of the KBE for finite systems |24[ I25| indicated 
unphysical long-time behavior (damping). Here, again solutions using single-time equations 
are, apparantly, closer to the exact result. It is, therefore, important to understand how the 
two-time and single-time approximations to the hierarchy are related to each other and whether 
and when single-time solutions are justified. 

The recipe how to derive the single-time approximation from the NEGF for an arbitrary 
selfenergy is the so-called generalized Kadanoff-Baym ansatz (GKBA) derived by Lipavsky, 
Spicka and Velicky |271 128] . The GKBA has, so far been used for spatially homogeneous systems, 
for electron-phonon scattering and for Coulomb scattering in the second Born approximation. 
Numerical comparisons for the case of Coulomb scattering have shown satisfactory agreement 
|29| 130] l3T] . Similar observations have been recently made for finite systems (26]. Nevertheless, 
a systematic analysis is still missing. The aim of this paper is two-fold. First, we discuss the 
relation between NEGF and DO (on the level of the second Born approximation). We take 
the opposite route, compared to Lipavsky et al. |27[ 128] and derive the GKBA result, starting 
from the BBGKY-hierarchy and identifying the approximations necessary to obtain the GKBA. 
Second, we consider, as an example of the quantum dynamics of a finite system, the evolution of 
a 4-site Hubbard nano-cluster — computed with Green's functions and the GKBA — and compare 
them to BBGKY-based results by Akbari et al. [JJ using the Wang/ Cassing approximation. 



2. Theory 

In this section, we give a brief overview of the theoretical foundations of the two aforementioned 
methods for the statistical description of quantum many-body systems. One of the most basic 
differences between these two approaches is the incorporation of the spin statistics of the particles, 
which, e.g., for Fermions leads to the effect of Pauli blocking |34J. While the Green's function 
method, which will be described in Sec. 2.2 has the spin statistics intrinsically build in by use of 
bosonic/fermionic creation- and annhilation-operators, the reduced density operator theory is 
formulated for spinless particles, which requires an explicit (anti-)symmetrization of the equations. 
This advantage of the former approach is at least partly counterfeited by the more complicated 
structure of the equations for the Green's function G(t,t'), which depends — apart from the 
physical time t — on an additional time argument t' , e.g., providing access to the corresponding 
ionized system. Numerically, this leads to comparatively more involved calculations with the 
Green's function approach. To make up for that, the introduction of the GKBA for the single- 
particle Green's function allows for a simple approximate reconstruction of the two-time Green's 
function from its time-diagonal value, by which a scaling of computation time comparable to 
density operator theory is achieved. 



2.1. Density operator theory 

For an iV-particle quantum system in a mixed state (in contact with a stationary environment) 
the proper fundamental quantity is the iV-particle density operator 
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where 



(fc) 

1...N 



is a possible micro-state compatible with a given macro-state of the whole system 



including the bath, and denotes the probability of its realization. This description is valid 
if the interaction of the system with the bath is weak. The density operator obeys the von 
Neumann equation 

ift%>i..jv- #i..jv , P1...N {i)=0, (2) 



denotes the standard 



where i?i...jv is the Hamiltonian of the TV-particle system and 

commutator. For a given initial state pi..jv(io) = Px*\n' * ne solution of this equation completely 
determines the time evolution of the system and is equivalent to the solution of the Schrodinger 
equation for all possible micro-states. 



2.1.1. Nonequilibrium quantum BBGKY -Hierarchy. Instead of the full iV-particle density 
operator, it is usually useful to consider reduced density operators, the s-particle operator F\,„ s 
being defined as a partial trace over the TV-body density operator p\...N P"T] . 



/•'....> = V^Tr,. :... v /':...v , ^Tr l ... s F 1 ... s = 1 



(3) 



where V denotes the volume, and the partial trace Tvi,,, s A of an operator A in an arbitrary basis 
of states \x) = \xi)\x2) ■ ■ ■ \%n) is defined as 
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The density operators obey a system of equations of motion — the BBGKY hierarchy that follows 
from taking the partial trace over Eq. ([2j 
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and so on. It is obvious that the whole hierarchy is equivalent to the von Neumann equation. 
Note that this system of equations is local in time, all functions depend on a single physical time 
t, and it has to be complemented by initial conditions Fi(to), -^12(^0) an d so on. Here Hi is the 
single-particle Hamiltonian whereas the two- particle Hamiltonian is given by H\2 = H\+ Hi + V^- 
As mentioned before, Eqs. ^ and ^ are written for spinless particles to make the mathematical 
structure more transparent. We will add the appropriate exchange contributions in the next 
section where we consider the (anti-)symmetrization and approximations to the BBGKY hierarchy 
via a cluster expansion. 



2.1.2. (Anti-)Symmetrization and Cluster expansion of the BBGKY hierarchy. For the practical 
analytical or numerical solution of the BBGKY hierarchy, it is obvious that it has to be truncated 
to become tractable. A suitable approach consists in performing the so-called cluster expansion 
to separate the two- particle, three-particle and higher correlations from the ideal part of the 
density operator which is given by a product of single-particle operators, 



F 12 {t) = Fx(t)F 2 (t) + ci 2 (t) , 

^123 (t) = F 1 (t)F 2 (t)F 3 (t) + F!(t)c2 3 (t) + F 2 (t)c 13 (t) + F 3 (t)c 12 (t) + c 123 (t) , 



(7) 



and analogously for the higher-order density operators. 

To correctly account for the spin statistics of bosons (fermions) we now ( ant i-) symmetrize 
all expressions. This can be done by introducing matrix representations with respect to an 
(anti-)symmetrized system of basis states defined in Fock space. Alternatively, (anti-)symmetric 
expectation values of observables can be computed with standard Hilbert space states when the 
respective operators are (anti-)symmetrized [35 . Here we follow the latter idea as it leads to 
more compact expressions. The (anti-)symmetrization of the density operators is achieved by 
replacing 
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where the binary/ternary (anti-)symmetrization operator is defined by its action on an arbitrary 
two-particle/three-particle state |12) and 1 123) , respectively, 
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and the upper/lower sign applies to bosons/fermions. 

We now introduce these cluster expansions ([7]) with the (anti-) symmetrized density operators 
into the BBGKY- hierarchy {[5],[6]) . We limit ourselves to the first three equations where the 
decoupling is achieved by neglecting four-particle correlations, C1234 = 0, (for details, see Ref. 
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Since below we will focus on the second Born approximation we do not explicitly write out 
the ladder and polarization terms |17| since they will be neglected. On the left-hand sides we 
introduced mean field Hamiltonians which are renormalized by a Hartree-Fock potential t/ HF , 

flj = H 1 + Uf F , Uf F = nTr 2 V 12 F 2 Af 2 , (14) 
H\2 = + H 2 , #123 = + H 2 + H 3 (15) 

and the non-hermitian operator V\ 2 = (1 ± nF\ ± nF 2 )Vi 2 which takes into account the exchange 
renormalization of the pair interaction giving rise, e.g., to Pauli blocking. 

Below we will consider two approximations: first, C123 = 0, which leads to the second Born 
(2B) approximation of NEGF theory together with the Hartree-Fock GKBA (HF-GKBA, i.e. 



with HF propagators), see Sec. 2.2 Second, we include in C123 all relevant terms that give rise to 
2B plus GKBA with full propagators (HF plus 2B). We start with the second approximation 
since the first follows from it as a special case. To this end, we solve the third hierarchy equation 
by retaining on the r.h.s. only terms that are proportional to c\ 2 : 

{ (^13 + ^23) F 3 =F nF 3 (F^rs + F 2 V 23 ) t n (c 13 V 13 + c 23 V 23 )} c 12 - h.c. (16) 

These terms describe the coupling of the pair 1-2 to third particles including medium effects 
which will give rise to energ renormalization (selfenergy) . This expression can be rewritten using 
the definition of V, and neglecting the ladder-type corrections involving products C13V13 and 
C23V23, yielding 

{(1 ± nFi) (1 ± nF 3 ) V 13 F 3 T F^Vn (1 ± nF 3 )} c 12 - h.c. + (1^2) 
=:(S? 3 TS? 3 + S>TS 2 %)c 12 -h.c., { j 

where we defined 

S%--=FiF*V ab F*. (18) 



We now rewrite Eq. (11) in static second Born approximation, neglecting all polarization and 



ladder terms and retaining, on the r.h.s., only the selfenergy contributions [17} 

±hd t C123 - {Hi 23 C123 - ci23#?23 fft } = ( S 13 =F S< 3 + S% 3 T S£ 3 ) c 12 - h.c. , (19) 

where we have introduced an effective (non-hermitian) three-particle Hamiltonian H^2 3 = 
Hi + H 2 + H 3 . We can now formally solve this equation for C123 in terms of C12, 

ci 23 (t) = U°+(tt )c 123 U^ 3 (tot) 

+ 4 f d ^ 2 + 3 (tt) { (^3 T ^3 + S 2 >3 T S 2 < 3 ) c 12 - h.c.} U° 123 (it), 

in Jtn T 
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where we introduced the propagators U® 23 with the properties 



U'Aitt') = [U$ 3 (t't)]\ U?f 3 (tt') = Ut{tH)Ut{tt')Uf{tt') . (21) 

The single-particle propagators obey effective one-particle Schrodinger-type equations 

\ihd t - fli(t)} U?{tt') = 0, Uf(tt)=Q, (22) 

ITf {t't) \ihd t . + H[(t') } = , C/f (fit) = , (23) 



where, in the second equation, the time derivative acts onto operators placed left of it. The 
definition of the effective one-particle Hamiltonian H\ will be derived in the following. 

We now turn to the second hierarchy equation (11), inserting the formal result for c\ 23 {t) on 
the r.h.s. and again neglect the polarization and ladder terms, 
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Here, the first term on the r.s.h. of Eq. (11) has been transformed according to 
V12F1F2 - FxF 2 V} 2 = (1 ± nFx) (1 ± nF 2 ) V 12 FxF 2 - F X F 2 V 12 (1 ± nF x ) (1 ± nF 2 ) 



where we introduced the greater- and less-collision integral operators, and If 2 
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Inspection of the integral term in Eq. (24) reveals that it has the structure of a selfenergy 
operator, £12, acting on c\ 2 : 



t 12 {t)c 12 {t) = J t dt {E+ 2 (tt)ci 2 (t)U^(tt) - h.c.} , 



which we can decompose into one- and two-particle contributions, 

S 12 = £1 + t 2 + tfj , ^ -^ ± TT± 1 ^±rr± , ^±,corr 



J 12 



EfUf + EfL/f + £ 



12 



(27) 



(28) 



which are, in turn, given by (cf. Eq. (24)), 
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With this, we can rewrite Eq. (24) by collecting all terms acting on c\ 2 into an effective 
two-particle hamiltonian, 



ihd tCl2 (t) - {H°f(t) Cl2 (t) - c 12 (t)H°f(t)} = im - Il2(t) , 

with the definition 
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This hamiltonian consists of three parts, H® 2 ^ = H\ + H 2 + Hy^ 01 , 
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Hi{t)c 12 {t) = fl?(t)cia(t) + f dtZf{tt)U 2 + {tt)c 12 {i)U° 2 {tt) , 

Jt 

H 2 (t)c 12 (t) = H°(t)c 12 (t) + / diU+(tt)E+(tt)c 12 (t)U^(it) , 



(32) 



HT(t)c l2 (t) = / dtX™ + (tt)c 12 (t)U» 2 (tt). 
J to 

To preserve the additivity of H? 2 resulting from the Born approximation, it is necessary to 
neglect the term H®2° r ci 2 , which yields a result that is consistent with our previous definition of 
the effective three-particle hamiltonian, 



H°fc 12 = (fli + H 2 ) c 12 , 

^123 ffc 123 = [H\ + H 2 + H 3 ) C123 



(33) 



This result can now be used to solve the equations of motion, Eq. ( 22 ) , for the renormalized 
one-particle propagator U^~(tt') which transforms into 



[ihd t - i?i j Uf(tt') - J* dt Et(l£)Uf(tt') = . 



(34) 



Using the definition of S% we can write out the selfenergy in second Born approximation 
explicitly (neglecting initial value terms), 



Sf(tt') = ^Tr 3 ^V 13 U+(tt')U+(tt') [F>F>V 13 F<tF<F<V 13 F>] U 3 (t't) 

The structure of suggests to define new quantities 

= U+(tt')F>(t') - F>(t)U-(tt') , 
9?tft) = U+(tt')F<(t') - F<(t)U-(tt') , 
giving rise to a compact and symmetric expresssion 
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where we defined the greater- and less-selfenergy Y^{tt') and Y,f(tt'). With these definitions one 
can finally write down the equation of motion for the single-particle density operator F± , which 
reads according to Eq. (10), 
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Here we have used, after the first line, the formal solution c\2{t) of Eq. (30) that includes initial 
correlations, c\ 2 and scattering contributions (the time integral). 

In the next section, we will see, that this equation of motion for the reduced single-particle 
density operator has the same form as the result obtained within the NEGF formalism for the 
time-diagonal limit of the function G < in second Born approximation, after applying the GKBA. 
Thereby it turns out that the full result, as derived above, corresponds to the GKBA with 
full second-Born propagators. In contrast, the HF-GKBA follows if the propagator equation 



(34) is solved after neglecting the selfenergy correction, which is equivalent to decoupling the 



BBGKY-hierarchy by requiring c\ 23 = 0, cf. Sec. 2.2 



2.2. Nonequilibrium Green's functions 

In contrast to density operator theory, presented in the preceding section, the basic quantities in 
the nonequilibrium Green's function formalism are the creation/annihilation operators at and 
o|, which create/annihilate a particle in the i-th. one-particle orbital and obey the canonical 
commutator /anticommutator relations for bosonic or fermionic particles, 
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With this, one can define the ensemble average of the combination N creation and annihilation 
operators at 2N different points in time, 
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as the less-part of the real-time iV-particle thermal Green's function G N,< . Here p is the density 
operator of the system and the notation "<" refers to the particular ordering of the operators. For 
equal times (t\ = . . . = tjs? = t' x = . . . = t' N ) the quantity i G J< is just the TV-particle density 
operator p\...N in the one-particle orbital basis. Analogous to the density operator theory, one 
can define reduced Green's functions involving fewer operators, so that the most basic quantity, 
the single- particle two-time less Green's function GfAt\,t2) can be defined as 



G<.(t 1 ,t 2 ) = (c J t (t 2 )Q(t 1 )) 



(41) 



The Green's functions also obey a hierarchy of equations of motion, the Martin-Schwinger 
hierarchy (MSH) [36J, which for all N connects the iV-particle Green's function to the (N — 1)- 
and the (N + l)-particle Green's functions. Since the complexity of the whole hierarchy — due to 
the additional time-arguments — is computationally even more demanding, one again resorts to 
the closure of the hierarchy equation on the single-particle level by a cluster expansion, stating 
the two-particle Green's function in terms of the one-particle Green's function by introduction of 
a suitable selfenergy functional S(G). With this, the equation of motion for the less Green's 
function attains the well known form of the Keldysh/Kadanoff-Baym equation (KBE), 

idt.Gf^h) = J2hik(tx)G<(t 1 ,t 2 )+ f dt 3 Y,f k (tx,t 3 )G<(t 3 ,t 2 )+ f dt 3 E&(ti,t3)G&(«3,t2). 
k J J 

(42) 

The quantities G R and G A — the retarded and advanced Green's functions — are defined as 

G* /A (h,t 2 ) = ±9(±(t 1 -t 2 )) [G^ihM) -Gf^hM)] , (43) 
G> ] {t u t 2 ) = -i(c i (t l )c){t 2 )) , (44) 



where hy is the ij-ih matrix element of the Hartree-Fock part of the Hamiltonian. The 
corresponding components of the selfenergy read in the second order Born approximation: 
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where a = {1,2} for spin-polarized/spin-restricted systems. As an alternative formulation of the 
KBE Lipavskii et al. proposed |27j 
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To further simplify this equation, they introduced the generalized Kadanoff-Baym ansatz (GKBA), 
which is equivalent to solving Eq. (46) in first order, 
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with the result 
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One notices that Eq. (47) is only formally closed in terms of the single- particle density matrix p, 
since the propagators G R ' A still obey a two-time equation of a similar degree of complexity as 
the original KBE. To circumvent this we introduce a further approximation, replacing the full 
propagators by the HF propagators GupS which are defined as 



G% A (t u t 2 ) = TiO[±(h - t 2 )\ exp \-i jf dt 3 h(t 3 ] 

Correspondingly, this approximation will be called HF-GKBA. 

Then the numerical solution of Eq. (47) can be obtained, using a finite time step A, in the 
following way |21j . 
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where 
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for a small time step A< 1, and the HF-Hamiltonian h is assumed not to change between t and 
t + A. The collision integral I(t, t) is given by 

J(M) = jf d * {s > '°(t,t)G < -°(t,t) -S < '°(t,t)G > '°(t,t)} , (53) 

where in all the two-time quantities under the integral the GKBA reconstruction according to 
Eq. (|48]) is used. 

When comparing Eq. ( 53 ) with the right-hand side of Eq. ( 38 ) , one immediately recognizes 
that it is of exactly the same structure, with the identifications 



E*'°(t,t) = E^(t,t), 

G R {t,i) = -e{t-i)u + {t,i) , 



(54) 



Moreover, the GKBA, Eq. (48) appeared naturally in our density operator theory, cf. Eq. (36). 
While this equivalence holds for general propagators U + defined by Eq. ( 34 ) , neglect of the 
renormalization (selfenergy term) — which is equivalent to decoupling the BBGKY-hierarchy 
by using C123 = — directly leads to HF-propagators and to the HF-GKBA of nonequilibrium 
Green's functions theory. 



3. Numerical example: 4-site Hubbard model 



To illustrate the findings obtained in sections 2.1 and 2.2 in this section we show some results 



for a 4-site Hubbard nano-cluster obtained from the Green's function method with the GKBA 
applied and compare to density operator results by Akbari et al. [1 . The Hubbard model 
|37} I3"8"] is a commonly used simplified description of a narrow-band solid-state system, where 
the motion of the electrons in the solid is mapped onto a hopping process between adjacent 
atomic sites with just one orbital for each spin projection. The Coulomb interaction between 
the electrons is assumed to be shielded so that it is mainly restricted to on-site interaction. 
Despite of these simplifications, this description is able to cover much of the rich behavior 
of these systems, for example the phase transition between a conductor and Mott insulator 
through the interplay of hopping and on-site interaction in two and three dimensions. Also in 
ID, the system is — depending on the interaction strength — strongly coupled and, therefore, the 
theoretical description needs to treat correlations between the electrons. 

The Hamiltonian of a one-dimensional Hubbard cluster comprised of N sites at electronic 
half-filling is given by 

N N 

H = ~tJ2J2 h V 4a £ ja + U 4t% g W ' ( 55 ) 
ij a i 

where c^J denotes the annihilation (creation) operator in the single-particle orbital on site i with 
spin a and hij = Suj\ is the hopping matrix between nearest neighbor sites with the convention, 
Snj) = 1, if are nearest neighbors, and <5{i,j) = 0, otherwise. The first term originates 
from the single- particle energies in the periodic lattice structure and incorporates the hopping 
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Figure 1. Density response on the first site n\ from strong nonequilibrium initial state for the 
4- site Hubbard chain with a coupling strength U = 0.1. The GKBA and the HF results are 
depicted by the solid red lines and green dashed lines, respectively. In blue the exact results from 
Akbari et al. [1] is shown. 



amplitude —t. The second term describes the on-site interaction of the electrons, which is given 
in terms of the interaction strength U. 

In the following, we consider a (N = 4)-site cluster with weak interaction, U = 0.1, at zero 
temperature with periodic boundary conditions. For this system we are interested in the study of 
strong non-equilibrium situations. To this end, we prepare the system in an initial state, where 
all the 4 electrons are forced to the left-most two sites and afterwards examine the free evolution 
of the system. For this setup, Green's function solutions within the HF-GKBA, as well as HF 
results are easily achieved. In Figs, [I] and [2j the time-evolution of the density n\ (T) on the 
first site — summed over the two orbitals for different spin orientations — measured in terms of 
the inverse hopping amplitude i _1 are shown. In Fig. [l] we present HF-GKBA results together 
with the exact solutions and time-dependent HF simulations. In Fig. [2] the HF-GKBA data are 
compared to exact results as well as to density operator results [1J within the Wang-Cassing 
(WC) decoupling of the BBGKY-hierarchy [TS]. 

The results show that the HF-GKBA performs very well for the propagation of a strong 
nonequilibrium initial state, especially in comparison with time-dependent Hartree-Fock, that 
does not show a decrease of the oscillation amplitude at all, present in the exact solution. 
Comparing the amplitudes with the exact ones, the GKBA overestimates them, though. In the 




Figure 2. Density response on the first site n\ from strong nonequilibrium initial state for the 
4-site Hubbard chain with a coupling strength U = 0.1. The GKBA results are depicted by the 
solid red lines. The blue pluses and brown dots represent the exact results and the ones obtained 
from the WC approximation to the BBGKY hierarchy by Akbari et al. pQ. 



part from T = 40 on, the GKBA cannot describe the oscillations sufficiently well, but always 
keeps a good phase agreement with the exact solution. Comparing with the WC approximation in 
Fig. [2] it is obvious that the WC solution has a very good agreement of the oscillation frequency 
up to T = 80, but the height of the peaks is much overestimated for time T > 40, both compared 
to the exact ones as well as to those computed with the GKBA. 

4. Discussion 

In this paper we have shown that, within the second order Born approximation, we can establish 
a one-to-one correspondence between nonequilibrium Green's functions within the GKBA and 
reduced density operator theory, yielding the same formulas for the equation of motion of the 
single-particle density matrix. On the example of the free evolution of a four-electron quantum 
dot from a strong nonequilibrium initial state we have shown the overall satisfactory agreement 
of the HF-GKBA to the exact dynamics of the system far away from equilibrium, for weak 
interaction strength. This is encouraging since the HF-GKBA allows for long propation times 
that are impossible to achieve with full two-time calculations and it retains most of the attractive 
properties of NEGF, such as time reversibility, total energy conservation and memory effects, e.g. 
|391I17|. In the future, it will be very interesting to examine the potential of the GKBA for higher 



order selfenergy schemes such as GW- or T-matrix (ladder) approximation. For the T-matrix 
approximation for homogeneous macroscopic systems, a similar correspondence has been found 
previously |4U[ II 7| . It will be interesting to extend this correspondence to inhomogeneous finite 
systems in order to achieve long propagations also for strongly coupled systems and to verify 
whether the GKBA is here able, as well, to remove the unphysical damping observed in two-time 
calculations [24]. 
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